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Abstract 

The electron and muon number violating muonium-antimuonium oscillation process in an 
extended Minimal Supersymmetric Standard Model is investigated. The Minimal Supersym- 
metric Standard Model is modified by the inclusion of three right-handed neutrino superfields. 
While the model allows the neutrino mass terms to mix among the different generations, the 
sneutrino and slepton mass terms have only intra-generation lepton number violation but not 
inter-generation lepton number mixing. So doing, the muonium-antimuonium conversion 
can then be used to constrain those model parameters which avoid further constraint from 
the /j. — > ey decay bounds. For a wide range of parameter values, the contributions to the 
muonium-antimuonium oscillation time scale are at least two orders of magnitude below the 
sensivity of current experiments. However, if the ratio of the two Higgs field VEVs, tan/?, is 
very small, there is a limited possibility that the contributions are large enough for the present 
experimental limit to provide an inequality relating tan/? with the light neutrino mass scale 
m v which is generated by see-saw mechanism. The resultant lower bound on tan f3 as a func- 
tion of m y is more stringent than the analogous bounds arising from the muon and electron 
anomalous magnetic moments as computed using this model. 

1 Introduction 

The time-dependent oscillation between two distinct levels or particle species is an interesting 
quantum mechanical phenomenon which has been widely studied in many physical systems vary- 
ing from a particle moving in a double-well potential of the ammonia molecule to oscillations 
in the neutral K° - K° and fi° - B° meson systems CD-El. It was suggested roughly 50 years 
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agoQ that there may be a spontaneous conversion between muonium and antimuonium resulting 
in an associated oscillation effect. Muonium (M) is the Coulombic bound state of an electron 
and an antimuon while antimuonium (M) is the Coulombic bound state of a positron and 

a muon (e + fi~). Since it has no hadronic constituents, muonium is an ideal place to test elec- 
troweak interactions. Of particular interest is that such a muonium-antimuonium oscillation is 
totally forbidden within the Standard Model because the process violates the individul electron 
and muon number conservation laws by two units. Hence, its observation will be a clear signal 
of physics beyond the Standard Model. Since the initial suggestion, experimental searches have 
been conducted! 5 ]-[ 6] and a variety of theoretical models have been proposed which can give 
rise to such a muonium-antimuonium conversion. These include interactions which can be medi- 
ated by (a) a doubly charged Higgs boson A ++ Q[8]], which is contained in a left-right symmetric 
model, (b) massive Majorana neutrinos l9l ITOl fTTTl . or (c) the r-sneutrino in an R-parity violation 
supersymmetric model lfT2l . 

In this paper we consider the muonium-antimuonium oscillation process in the Minimal Super- 
symmetric Standard Model extended by the inclusion of three right-handed neutrino superfields. 
While the neutrino mass terms can mix inter-generationally, we allow only intra-generation lep- 
ton number violation but not inter-generation lepton number mixing for the sneutrino and slepton 
mass terms. In this model, there are intermediate states which can contribute to the muonium- 
antimuonium oscillation process but not to the /u — > ey decay. Therefore, the experimental limit 
on muonium- antimuonium oscillations can be used to constrain those model parameters which 
are not constrained by the /j, — > ey decay bounds. In order for there to be a nontrivial mixing 
between the muonium and antimuonium, the individual electron and muon number conservation 
must be violated by two units. Such a situation will result provided that the neutrinos are massive 
Majorana particles or the mass diagonal sneutrinos are lepton number violating scalar particles. 

2 The extended Minimal Supersymmetric Standard Model 

The Minimal Supersymmetric Standard Model (MSSM) is the supersymmetric extension of the (2 
scalar doublet) Standard Model with the minimal particle content! 13]. For each particle, there is 
a superpartner with the same internal quantum numbers, but with spin that differs by half a unit. 
Tab. 1 lists all the chiral supermultiplets needed for MSSM, 
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Names 


Spin 


Spin^ 


SU(3) c ,SU(2) L ,U(l)y 


squarks 
and quarks 
(x3 families) 


Q 

u c 

D c 


(u L d L ) 

— 

u R 

3* 

R 


("La dux) 
R 

d" 

R 


(3,2,1) 
(3,1,-f) 
(3, 1, 1) 


sleptons, leptons 
(x3 families) 


L 

E c 


(vl e L ) 

e* 

e R 


(VLa eLa) 
4 


(1,2,-1) 
(1,1,1) 


Higgs, higgsinos 


Hj 
H B 


(h+ h° T ) 


(Ka ~ h rJ 


(1,2,1) 
(1,2,-1) 



Table 1: Chiral supermultiplets of the MSSM 



while Tab. 2 summarizes the gauge supermutiplets of MSSM. 



Names 


Spin \ 


Spin 1 


SU(3) c ,SU(2) L ,U(l)r 


Gluino, gluon 


la 




(8,1,0) 


winos, W bosons 


"a "a 




(1, 3, 0) 


bino, B boson 


B a 




(1,1,0) 



Table 2: Gauge supermultiplets in the MSSM 



In the above two tables, the dotted and undotted indices, a, a, indicate 2-component Weyl spinor 

fields. In the subsequent analysis, we will recast all the spin 1 fields as 4-component Dirac spinor 

fields, which will be represented using the same symbols, but without the dotted and undotted "a"s. 

( \ 

v La 

For example, VLa is the Weyl representation of the left-handed neutrino field, while vl = 
is the Dirac field. 

In order to implement the see-saw mechanism[14] for neutrino masses, we consider an exten- 
sion of the MSSM, where one adds three additional gauge singlet chiral superfields N? (i=e,/u,r 
denotes the generation), whose ^-component is a right-handed neutrino field, 

**i = *«O0 + ^<?*VR(y)a + (POaFwiy), (1) 

where 

/ = ^ + ie*a^S*. (2) 
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These S U(3) x SU(2) L x U(l) singlet superfields are coupled to other MSSM superfields via the 
superpotential. We employ the most general R-parity conserving renormalizable superpotential so 
that the superpotential is 

W = -iue ah H a B H h T + AieabEfflH* + X^H^N) + ^M^N c p (3) 

while the relevant soft supersymmetry breaking terms are 

rEMSSM t 'j\2/~* ~ , ~p* ~h \ i U\2p* t < >j\2~* ~ 

£ S0 ft = -t m D \ v iL v jL + tutjL) ~ (mp { iR tjR - (m$ v iR v jR 

-(XijAij^ViLvfr + MpijVmVjR + AiCuhltaJt^ + C«X4/* R + H.C.). (4) 

The interaction terms that contribute to the muonium-antimuonium oscillation and the electron 
and muon anomalous magnetic moments can be extracted from the Lagrangian of this extended 
Minimal Supersymmetric Standard Model (EMSSM) as 

JZ = —(w-^iLrMViL + W + »v iL y^ iL ), (5) 

Ll t = -igi^iL W~ v iL -V iL W £a), (6) 
^ = ^iL^l iL -l* L Wn iL ), (7) 

£f nt = ^(iu. Bl iL -7* L B £ iL ) + yj2 gl i(f iR B~l iR -7* R B f iR ), (8) 
_ ij t £/\* 

4 = <« + *J fc) - ^ iR h° B 7 iL + 7* L W B e iR ). (io) 

In the above equations, all the spin ^ fields are Dirac spinor fields. In particular, note that the field 
h B has the Weyl field decomposition 



h B = 



( T- \ 
h + 



(11) 



The parameters V R and Vj are the vacuum expectation values of the two Higgs fields: < h B >= Vb 
and < hj >= Vj. These VEVs are related to the known mass of the W boson and the electroweak 
gauge couplings as 

2M 2 

Vl + V$ = V 2 = — ^ * (\14GeV) 2 , (12) 
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while the ratio of the VEVs is traditionally written as 



tan/3 = — . 
V B 



(13) 



In the above, - X'^Vt is the Dirac mass matrix of neutrinos and m, are the lepton masses. 
Since the masses of electron and muon are small, the terms which have couplings proportional to 
nti/Vs in interactions (9) and (10) are severely suppressed and will be ignored in the subsequent 
analysis. 

The neutrino mass term can be extracted from the superpotential terms A'^abH^L^Nj and 
jM i [ >NfN c j as 



■Cmass = ~\( ( V L) C Vr ) 



m 

m D 



D 

Mr ) 



vl 
(vrY ) 



+ H.C., 



where 



vl = 



VeL 




/ \ 

VeR 




vr = 


V^R 


. V rL , 




. V rR , 



(14) 



(15) 



Note that the upper left 3x3 block in the neutrino mass matrix is zero. This block matrix involves 
only left-handed neutrinos and in our EMSSM its generation requires a nonrenormalizable super- 
potential term. Consequently we ignore this term. For three generations of neutrinos, the six mass 
eigenvalues, m va , are obtained from the diagonalization of the 6 x 6 matrix 



M v 







m, 



m D M R 

Since M v is symmetric, it can be diagonalized by a single unitary 6x6 matrix, V, as 

Kia g = V T M V V. 

This diagonalization is implemented via the basis change as following 



(16) 



(17) 



vl 
(vrY 



l \ 

VeL 




( \ 
V\ 


V M L 




v 2 


V T L 


= V 


V3 






{VeRf 




(v 4 ) c 


(v m rY 




(v 5 ) c 









vr 



' (v eL y ' 




' (viY ' 


(v ML ) c 




(v 2 Y 


(VtlY 


_ y* 


(V 3 ) C 






VeR 




y 4 


V M R 




V5 


, V tR , 




. y 6 , 



(18) 
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The neutrino mass term then takes the form 

1 6 _ — 



m va v a v a , 



(19) 



0=1 



a=l 



where m va are the Majorana neutrino masses. 

Since a nonzero Majorana mass matrix M 1 ^ does not require SU(2)l X £/(l) symmetry break- 
ing, it's naturally characterized by a much larger scale, Mr, than mo, which is the scale of the 
Dirac mass matrix m 1 ^ whose nontrivial value does require SU{2)i x U(l) symmetry breaking. So 
doing, one finds on diagonalization of the 6x6 neutrino mass matrix that the three eigenvalues are 
crudely given by 



111, 



«: inn, a = 1,2,3, 

Mr 



while the other three eigenvalues are roughtly 



Illy 



Mr, 



4,5,6. 



(20) 



(21) 



This constitutes the so called see-saw mechanism lfT4l and provides a natural explanation of the 
smallness of the three light neutrino masses. Moreover, the elements of the mixing matrix are 
characterized by an mo /Mr dependence 



Vab-Oil), a,b = 1,2,3, 
Vab~0{\), a, b = 4, 5, 6, 

a = 1,2,3,6 = 4,5,6. 



„ mn 
Vab ~ V ba ~ O , 
Mr 



(22) 



In order to obtain the sneutrino masses, it's convenient to define vn - -^=(v ( xi + ivm) and 
Vir = -^(v,t?i + ivjR2). Then, the sneutrino-squared mass matrix separates into CP-even and CP- 
odd blocksIBl, 



■Cmass - 2 ( ^1 ^2 ) ^ 



i,J=e,[i,T 



( <t>i 



2 / 





4} 



(23) 



where <f>' a = {vn a ViR a ) and Ai?. consist of the following 2x2 blocks: 



(m ! /) 2 + ±m§ cos 2/3 + (m^) 2 



m^(A, 7 -//cot^±M^) 



m^CAy - H cot/? ± M l R J ) (M^) 2 + (m^) 2 + (m^) 2 ± 2%M£ 



(24) 



6 



with Ajj and Bij are SUSY breaking parameters (cf. Eq.(4)). Since we allow only intra-generation 
lepton number violation but not inter-generation lepton number mixing for the supersymmetric 
partners, we can arrange the parameters in matrices (24) so that At? , + = for i + j. So doing, the 
sneutrino mass term simplifies into three 4x4 matrices for three generations 



■Onus = J] 1( 4 ) 



i=e,fi,r 



Ml 







Mi 



<t>\ 



(25) 



2 J 



The sneutrino mass matrix M\.. is real and symmetric, so it can be diagonalized by a real orthog- 
onal 4x4 matrix, U l , as 



M 



Vjdiag 



= U iT MlU\ 



where U' is in a form as 



U l = 



' u\ o 1 

Ul 



(26) 



(27) 



This diagonalization is implemented via the basis change on <p\ and (f>' 2 



where y,y, are all real. Then the sneutrino mass term takes the form 





( ~ \ 

ViLl 




( ~ \ 

Vfl 


\ 




= U l 


va 


/ 


ViL2 




Vi3 




K ViR2 , 







(28) 



i 4 

r vi _ _£ V ^~ ~ 



(29) 



a=\ 



where m y J are the sneutrino mass eigenvalues. 

In the following derivation we assume that M n R is the largest mass parameter. Then, to the first 
order in 1 /M R , the two light mass eigenvalues are roughly 



in *K) 2 + > Z COs2^- 



1 



2(mil) 2 (A i7 -ncotfi-Bu) 



m~ « (nil) 1 + cos 2/3 + 



1 



limlfiAu-LicotB-Bu) 



M" 



while the two heavy mass eigenvalues are 



M R 



(30) 



7 



"4 * W - 2B " m r- 



f ii\2 



(31) 



To avoid excessive complication in our calculations, we expand U' in powers of the matrix 

trii 

parameter ft = The form of U to first order of ft is 



C/«. 

o uL 

i ft 
-ft i J 







I ft 1 J 



(32) 



The slepton mass term is given by 

■C-mass ~ ( ^iL ^iR ) 



i,j=e,fi,T 



{mf) 2 {mf) 2 



1 t ^ 



where 



{mf) 2 = ( m ! /) 2 + ml cos 2/?( sin 2 8 W - -), 
{mf) 2 = A if iV T + AidiVB, {mff = A^Vb for i * j, 



{mf) 2 = (ml) 1 - mi cos 2/3 sin z 6 W . 



(33) 

(34) 
(35) 
(36) 



In analogy to the sneutrino mass term, we can arrange the parameters in Eq. (34)-Eq.(36) so 
that (m^) 2 , {mf) 2 and {mf) 2 are zero for i + j and there is no inter-generation lepton number 



c, 



mixing in the slepton mass term. So doing, the slepton mass matrix reduces to three individual 
mass matrices for three generations 



■E-mass ^ ( ^ir ) 



i=e,fi,r 



{mf) 2 {mf) 



LR\2 



{m^) 2 {m™) 1 I 



RR\2 



(37) 



Since /l,-Vg = m,-, the off diagonal matrix element, {mf) 2 , can be written as 



{mf) 2 =m i {n\anp + C ii ). 



(38) 



Because the masses of electron and muon are very small compared with the sparticle mass scale, 
we ignore these off diagonal terms and consider 4l and (m as mass eigenstates. 
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Inserting the transformation (18) and (28) in the interaction terms (5)-(9) yields the explicit 
interactions in their mass basis: 



' i=e,fi,T a=l 



(39) 



-C = -^ Z Z^'^+^l Z Zfow-i/Lv fa +».c., (40) 



V2 



i=e,/x,T a=l 



V2 Z 

* i=e,fi,T a= 



i=e,n,T 



2 y 
mi 



i=e,fi,T 



4 • 17 

im T 



(41) 



4L = ^ Z B ^ - ^ B ^) + ^ Z ('* B - ^ B '4 (42) 

' i=e,/i,T 



£ int= Z Z^^^^ Z Z^-^*i^ +aC " (43) 



3 The muonium-antimuonium oscillation in the EMSSM 

The lowest order Feynman diagrams accounting for muonium and antimuonium mixing are dis- 
played in Fig.l. Graphs (a), (b), (c) and (d) are the non-SUSY contributions, which are mediated 
by Majorana neutrinos and W boson. The other graphs all involve SUSY partners. Graphs (e) and 
(f) are mediated by sneutrinos and winos, while graph (g) and (h) are mediated by sneutrinos and 
higgsinos. 





(a) 



(b) 
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Figure 1: Feynman graphs contributing to the muonium-antimuonium mixing. 



The T-matrix elements of graphs (a) and (b) are lfTTTl 



4 



T a = T b = -—J-— m 3) 7tl (l - r5 )e(2)][/i(4)y(l - y 5 )e(l)] 



256tt 2 M2, 

6 



(44) 



• #=1 a,b=l\a±b 

where /i(3) = //(/?3> ^3) , /i(4) = ju(p4, 54) , ^(1) = e(p\, s\) and e(2) = e(p2, Sz) are the spinors of 

2 

the muons and electrons and x Ya = , a = 1, 2, 3, ...6. We define S(x Vo ) and r(x Va , x Vi) ) as 

x 3 -llx 2 + 4x 3x 3 
4(1 - x) z 2(1 - xy 



,R(x A ) - R(x B ) 
1 (X A , x B ) = x A x B \ 

with 



T(x A , x B ) = x A x B ( — — — ) = T(x B , x a ), (46) 

v x A -x B ' 



R (x) = x2 - Sx+ ? 4 lHx)- 3 -^—. (47) 
V ' 4(1 - x) 2 K ' 4 (1 - x) 



The T-matrix elements of graphs (c) and (d) are lfTll 

T c = T d = 256n l M 2 ^(3)^(1 - 7s)«(2)] [^(4)^(1 - 7s)e(l)] 



{£(^a0 2 ^v )+ ^ (^ 2 (^) 2 tf(x Vfl ,x Vi )]. (48) 

a=l ai>=\;a+b 



The functions G(x v ) and ^(x v , x v , ) take the forms 



with 



(x A -4)x A (x 3 -3x 2 +4x A +4)x A 

(x A - iy 2(x A - \y 



v , x , L(x A ,x B ) - L(x B ,x A ) 

K(x A ,x B ) = yJx A x B , (50) 

x A - x B 



Tl N 4 - x A x B x A (2x B - x A x B - 4) 

L(x A ,x B ) = — — + — lnx A . (51) 

2(x A - 1) 2(x A - l) 2 



The T-matrix elements of graphs (e) and (f) are 

T e = Tf= - r^TTTTo [fiO)7 M (l ~ 7sM2)]&i(4)7p(l - ys)e(l)] 



10247r 2 M?- 
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/ 2 2 2 4 



£1=1 6=1 £1=1 fc=3 

4 2 4 4 , 

Z Z^^'OV**) + Z Z^sa)^^) 2 ^,^) , (52) 



where 



£i=3 b=l £1=3 fo=3 



? 



(53) 



with 



/(xi) - 7(x 2 ) 

/(*i,*2) = , (54) 

x\ - x 2 

x 2 In x - x + 1 

J M = — 777— ■ (55) 

(x - l) 2 



Finally, the T-matrix elements of graphs (g) and (h) are 

Z(n^niT){m e J.m e J.*) 

i,j=e,n,T T h - 

2 2 ' 2 4 



Z Z (C/ 2£i) 2 (^) 2/ fe !fl '^,) - z Z (f/ ^ )2(f/ i )2/fe -^,) 

£1=1 i)=l £1=1 fo=3 

4 2 4 4 \ 

Z Z ( ^£i) 2 (^) 2/ fe^) + Z TF\a)\<f^ la ,Z,A (56) 



£1=3 b=\ a=3 b=3 



where 



mi. 

Zy ia = -77T. (57) 
Mi 



4 The effective Lagrangian 

Combining all the T-matrix elements, we secure an effective Lagrangian which can be cast as: 



Leff = -ZZ-lftfV ~ 75)e][fi7^l - y 5 )e], (58) 
V2 



where 



4 / 6 



V2 5\2n 2 M 2 



W X £1=1 



Z(^0 2 (5(x y J-G(x v j) 



+ Z ((v,av: a )(v Mh v: h )T(x Va , Xvb ) - (v, a f(v: h ) 2 K(x Va , Xvb ))\ 

a,b=\;a+b ' 
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4 / 2 2 2 4 

52 



- 2048 J M 2 • 2 2 ( ^ )2( ^ )2/( ^-^ ) - 2 Z^^ 7 ^.^) 

jy- ' a=l Z>=1 <3=1 b=3 

4 2 4 4 v 

- 2 2<^ )2 w^wo + 2 2<^ )2 w 2/ ^'M 

a=3 6=1 a=3 6=3 ' 

y K< ) / y y & 2 j 2 

i,j=e,/x,T T h - a =\ b=\ 

2 4 " 4 2 

- 2 2 (£/ 2 fl )2(£/ i )2/(z ^' z ^ ) - 2 2 (C/ 4«) 2 (^) 2/ fe, a ,zv ; ,) 

a=l 6=3 £1=3 2>=1 

4 4 \ 

+ 22 (C/ 4£ 1 ) 2 (^i) 2/ fe M '^) • (59) 

£i=3 b=3 I 

Muonium (antimuonium) is a nonrelativistic Coulombic bound state of an electron and an 
anti-muon (positron and muon). The nontrivial mixing between the muonium ( \M > ) and an- 
timuonium (|M >) states is encapsulated in the effective Lagrangian of Eq. (58) and leads to the 
mass diagonal states given by the linear combinations 

\M ± >= 1 [(1 + s)\M > ±(1 - e)\M >], (60) 
V2(l + M 2 ) 

where 



£ = , : , (61) 

V^MM + V %M 

<M|-p 3 r£ e// |M> <M|-p 3 r£ g// |M> 
M M m - — ; _ _ Mmm - — ; _ _ (62) 

V< M\M >< M\M > V< M\M >< M\M > 

Since the neutrino sector is expected, in general, to be CP violating, these will be independent, 

complex matrix elements. If the neutrino sector conserves CP, with \M > and \M > CP conjugate 

states, then M MM = M MM and e - 0. In general, the magnitude of the mass splitting between the 

two mass eigenstates is 



\AM\ = 2 



ReyjM M MM 



MM 



(63) 



Since muonium and antimuonium are linear combinations of the mass diagonal states, an initially 
prepared muonium or antimuonium state will undergo oscillations into one another as a function 
of time. The muonium-antimuonium oscillation time scale, t^ m , is given by 

— !— = |AAf|. (64) 

T MM 
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We would like to evaluate \AM\ in the nonrelativistic limit. A nonrelativistic reduction of the 
effective Lagrangian of Eq. (T58l) produces the local, complex effective potential 

Veff (r) = 8%^ 3 (r). (65) 
Taking the muonium (anitmuonium) to be in their respective Coulombic ground states, 0ioo(r) = 

i ~ i in in 

e~ r/a MM^ where a^ M = — ! — is the muonium Bohr radius with m rec i - — — - m e the re- 



4 



. 3 ' tyijyi m da 

MM 



duced mass of muonium, it follows that 
1 



J d 3 rf m (r)\ReV ef f(rMr) 10 o 



T MM 

\ReGgM\ . , fl 16 \ReGMM\ 1 

16 = — |0ioo(O)| - — 5 — ■ (66) 

V2 n V2 a 3 . 



MM 



Thus we secure an oscillation time scale 



1 ^6\ReG MM \ my (67) 



T MM n V2 

5 Estimate of the effective coupling constant 

The present experimental limit] 6] on the non-observation of muonium-antimuonium oscillation 
translates into the bound 

\ReG MM \ < 3.0 x KT 3 G F , (68) 

where Gp =; 1.16 x \0~ 5 GeV~ 2 is the Fermi scale. This limit can then be used to construct some 
constraints on the parameters of this model. 

For simplicity, we set the neutrino Dirac mass matrix elements ni 3 D and the right-handed neu- 
trino mass matrix elements M*i to some common mass scales mp and Mr respectively. The light 
neutrino mass scale m v is of order m 2 D /Mn, while the heavy neutrino mass scale is of order Mr. 

Using these assumptions and taking into account the mixing matrices approximations Eq.(22) 
and (32), we can simplify the effective coupling constant Eq.(59) to a more manageable ap- 
proximated form. The contribution from graphs (a), (b), (c) and (d) in G^ M is ~ 5n Jt M 2 ' 

W 

(2ti(v^C) 2 (s(*vJ-G(*vJ^ 

With the limits of m V] , m V2 , m V3 ~ 0(jf-) and m V4 , m Vs , m V6 ~ 0{Mr), the contribution of graphs 
(a), (b), (c) and (d) can be approximated as 

easel: — — - In — , a=l, 2, 3, b=l, 2, 3, 

5nn 2 M 2 w M 2 R M 2 W { m 2 D ) 
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case 2: — • In 5 — , a=l,2, 3, b=4, 5, 6, 

5\2n 2 Ml M 4 R M 4 W { m 2 D ) 

case 3: • -^r- In , a=4, 5, 6, b=4, 5, 6. (69) 



SYln 2 M\ M\M 2 W \M W 

Taking Mr as the largest mass parameter, the first case and the third case are comparable, while 
the second one is suppressed by a factor m^/(M|M^). Therefore, the contribution from graphs 
(a), (b), (c) and (d) is roughly 



st 



Y,{V, a V: ci ) 2 (S{x Va )-G{x Va )) 



5\2n 2 M 2 , 

W V £1=1 

6 

+ 2 {( v yaV; a )(VKbV: b )T(x Va ,x Vb ) - (V Ma ) 2 (V: h ) 2 K(x Va ,x Vb )) 

a,b=\\a+b 

256;r 2 M|M^ M w 

The second term in Eq.(59) is the contribution of graph (e) and (f), in which the function I(y VlJa ,yv eb ) 
is a decreasing function of y$ and y Veb . It will be small for heavy seutrinos. To see this, we employ 
the approximations Eq.(32) 



u tv u \v u iv U^-CKl), 

U\ v U{ 2 , U^-Oi^X (71) 



so that the terms involving heavy sneutrinos will get an extra suppression from the mixing matrix. 
Therefore, the contribution of graph (e) and (f) is dominated by the term that only includes the 
light sneutrinos so that 

- 20 jj M 2 ■ (Z IX) W'0V'3o - 2 ix^^'ov^) 

yj- ' a=l b=\ a=l b=3 

4 2 4 4 

a=3 b=l a=3 b=3 

#2 

~ ~ 204Sn 2 Mi ' v 7 ^" 3 ^ ~ - ^(yv p3 .yv,i) + ^(yv^.y^))- ( 72 ) 

w- 

Employing the squared-mass difference between the two light sneutrinos in Eq.(30), the above 
expression can be approximated as 

( 7 0%, . yvei ) - Ky v ^ ) - 7 0%3 . Jv el ) + / 0% 3 , yv ei )) 



204Sn 2 M 2 - 
w 



15 



4 . _ v ,_d d 

4 Am j A *4 _i c* 

* 2048^M 2 _ 'Ml ' Ml ' dy^dyj^^ Ui) 

w- w- w- ^ el 

where the squared-mass differences are 

4(m^) 2 (A w -//cot^-B w ) 
, 4(m" ) 2 {A ee - u cot /? - 

K - J. — -• (74) 

Assuming A m - A ee = A and B m - B ee = B, the squared-mass differences of light muon 
sneutrinos and light electron sneutrinos are 

„2, 



4m£ ) (A- J ucot/?-B) 
Am?, - Am 2 £ = Am 2 = — 



(75) 



so that Eq.(73) then simplifies to 



g\ Am i Am 2 d q g^mj) 2 d d 



" 2048^M|. ' M|. ' M|_ ' ^ ft*/ *^ * 2048^M|. dy^ dyj^^ (76) 

The contribution from graph (g) and (h) is not dominated by the terms involving only light 
sneutrinos even though I{Zv ia ,z Vjb ) is a decreasing function of z Via and z Vjh , because these terms get 
suppressed by the mixing matrix. The terms including only light sneutrinos v,i, Vj\, vp 

roughtly gives 



(U i 43 ) 2 (U J 21 ) 2 I(Zy ,Zy fl ) + {U^f(U^) 2 I{Z Va ,Z Vfl ) 



(U l 2l ) 2 (U J 2[ ) 2 I(Zy n ,Z Vjl ) - (U 2l ) 2 (U J 43 ) 2 I(Z Vn ,Zy j3 ) 

:u\j 2 {u j 2i ) 2 i{z VB ,zy fl ) 
J> d d 

" Mg M\_ dz, a dz Vj r Z ^ l) 

while the terms including one light and one heavy sneutrino are roughly 

(U l 2l ) 2 (U J 22 ) 2 I(zv n ,z Vj2 ) - (^ 1 ) 2 (t/i 4 ) 2 /(Zv,,Zv, 4 ) 
AU^)\U j 22 ) 2 I(zy a ,z Vj2 ) + {U^) 2 {U{ A fl{z VB ,z V] ,) 
+(U' 22 ) 2 (U J 21 ) 2 I(z Vl2 ,z Vjl ) - (U l 22 ) 2 (U J 43 fl(z Vl2 ,z Vj3 ) 



(77) 
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m D 
Mr 



2 AM? - Am? (Mi-^ A 



m D 



M~ 



-3-)". 



2 AM? - Am? /-M,-- 



M~ 



Mr 







' 1 ^ 



M 6 



(78) 



where AM? is the heavy sneutrino squared-mass difference 

AM? = 4B u Mi 

Under our approximations, 



(79) 



AM? = AM? = ABMr. 



(80) 



The terms including two heavy sneutrinos are roughly 



(U l 22 ) 2 (U J 22 ) 2 I(z Vl2 ,z Vj2 ) - (U 22 ) 2 (U J 44 ) 2 I(z Vl2 ,z Vj4 ) 
-(U 44 ) 2 (U J 22 ) 2 I(z Vl4 ,z V]2 ) + (U 44 ) 2 (U J 44 ) 2 I(z Vl4 ,z V]4 ) 
AM? ■ AM? M\_ 



M- 
h B 

(AM?) 2 M? 

h B 



3M% 



3M| 



(81) 



Comparing the Mr dependences of Eq.(77), (78) and (81), we see that the dominant term is the one 
involving two heavy sneutrinos. Thus the contribution from graph (e) and (f) can be approximated 

as 



i,j=e,fi,T 
2 4 

V Wrr« \2(ttJ ^2 



2 2 

V V/rr« \2,r r ./ ^2 J 



2048VW \hU Ula) iU2h) I(Z ^ b) 



£1=1 fe=l 



4 2 

2, Z (f/ 2a) Z (^) Z/ fe ia ,^,) - ^ 2 (C/ 4 fl ) 2 (^) 2/ fe, a '^) 

£1=1 Z>=3 £1=3 b=l 

4 4 

.2 (TJ i x2, 



a=3 fo=3 

3m* (AM?) 2 3g^m4(AM2) 2 (l +tan 2 /3) 2 



2048V^2 M 6 



8 192/r 2 M|M^, tan 4 /? 



(82) 
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Combining the various contributions, the effective coupling constant is thus roughly given by 
\ReGft M \ « 



' ' -in — - — • - - Kyp^yfei) 



2WM\M% M w 2048/r 2 Mt_ dy^dy^ 
34m4(AM|) 2 (l+tan 2 ^) 2 

(83) 



8l92n 2 M 6 R M^ tan 4 p 

The first term in Eq.(83) is the dominant contribution of graph (a), (b), (c) and (d), which contains 
the intermediate neutrino and W boson. This contribution appears in the model in which all SUSY 
partners decoupled. The second term is the dominant contribution of graph (e) and (f), in which 
wino and sneutrino appear in the intermediate states. Finally, the third term is the dominant contri- 
bution of graph (g) and (h), with intermediate higgsino and sneutrinos lines. The second and third 
terms both depend on the sneutrino mass splitting. This reflects the intra-generation lepton number 
violating property of the muonium-antimuonium oscillation process, because the sneutrino mass 
splitting is generated by the AL = 2 operators in the sneutrino mass matrix. To compare the rela- 
tive sizes of these three terms, we use the current experimental limits of the neutrino and sparticle 
masses. 

The first terms in Eq.(83) has a factor m 4 D /M 2 , which is the scale of the light neutrino mass 
square m 2 - m 4 D /M R generated by see-saw mechanism. The experimental constraints on neutrino 
masses are summerized in reference |[T9l as 

m v (electron based) < 2eV, 
m v (muon based) < 0A9MeV, 

m v (tau based) < \S.2MeV, (84) 



For instance, assuming 



m D ~ M w , (85) 



2 

mi 



m v = ~ leV, (86) 

M R 

then the right-handed neutrino mass scale is about 

Mr ~ \Q n GeV. (87) 

In this case, the first terms in Eq.(83) is roughly 

9gim 4 n Mr 9gjml Mr 9R , 

2 f -. ■ In — = — ^-V-hi— » 3.9 x lO- 28 GeV~ 2 . (88) 

256n 2 M 2 R M A w M w 256n 2 M A w M w 
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The second term in Eq.(83) depends on the light sneutrino squared-mass difference Am?, which 



can be written in terms of light sneutrino mass splitting Am v by 



Am v = 2m v Am v , 



(89) 



where m v is the mass scale of light sneutrinos. So doing the second term in Eq.(83) can be written 

as 

g 4 2 (Amj) 2 d d d d 

-Hyv^yvei) = - — TTTf. ■ ~ ~ i Cy^ui > yv e i ) ■ (90) 



204Sjt 2 M 6 - dy v . dy 



'ill "JVcl 



5\2n 2 M 6 - dy- v . dy 



w- w- 
Reference IPT51 provides an upper limit on the sneutrino mass splitting by calculating the one-loop 
correction to the neutrino mass. Assuming that this correction is no larger than the tree result gives 



Am v < 2 x 10 3 ra v . 



(91) 



Relaxing this absence of fine tuning constraint can substantially enhance the contribution of the 
graph (e) and (f). Taking the sneutrino mass splitting to be of the same order as sneutrino mass 



Artiy ~ nvf, 

and rriy , m Ve to be the common mass scale m v gives 



2 



(92) 



(93) 



Eq.(90) can then be written as 

g\m 2 -{Am v ) 2 _d d_ 

~ 512jt 2 M^ ' dy Vtil dy % , 

4 

The function -jjg- ■ * ^ I(y Vul , y %l ) 



4 4 
Si m v 



d d 



5\2n 2 M 6 - dy v . dy 



w- 



ixl "JVel 



(94) 



'ill "J y e\ 



in Eq.(94) is a decreasing function of My,_. In 



order to calculate the maximum contribution of graph (e) and (f), we use the experimental lower 
bound on M^,_ . Many experimental searches for physics beyond the standard model have been 
conducted and provide various constraints on SUSY parameter space. Tab.3 lists some of the 
constraints! 18]. 



Sparticle lower limit [GeV] 


Sparticle lower limit [GeV] 


ft 94.0 
*? 46.0 
v 94.0 


fin 94.0 
e L 107.0 
e R 73.0 
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Table 3: Experimental lower limits on SUSY particle masses 



Fixing the wino mass to its lower limit in Tab. 3 

My- = 94.0GeV, (95) 



4 4 

the contribution • ^-^-I(y^,y,J 

is shown in Fig.2. 



as a function of sneutrino mass scale m v 




Figure 2: The contribution of graph (e) and (f) as a function of m v when fixing the wino mass to 
its lower limit M^>_ . 



When niy = 94.0GeV, which is allowed by the experimental limit in Tab.3, the contribution of 
graph (e) and (f) reaches its maximum so that 



*X d d 



< 1.3 x \Q- l0 GeV- 2 . (96) 



Finally, the third term in Eq.(83) depends on the heavy sneutrino squared-mass difference 
AM? = 2MyAMy = ABM R . Since we assume that Mr is the largest mass scale, AMy can't be 
arbitrarily large. Taking parameter B one order of magnitude smaller than Mr, the heavy sneutrino 
mass splitting is 

AM, ~ ^. (97) 
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The contribution of graph (g) and (h) can then be written as 

3g 4 m 4 (AM 2 ) 2 (1 + tan 2 /?) 2 3g 4 m 2 (l + tan 2 ^) 2 



(98) 



S\92jr 2 M 6 R M^tan 4 /3 204800/^^ tan 4 /?' 

When tan/? is very small the contribution can get large and even reach the experimental limit 
Eq.(68). In this case, the experimental limit of muonium-antimuonium oscillation provides an 
inequality relating tan/? and m v , which is given by 



3s 4 m 2 (l + tan 2 /?) 2 



< 3.5 x \0-*GeV- 2 . 



2048007r 2 M 4 v tan 4 J S 

This inequality translates into a lower bound of tan/? for different light neutrino masses m v \ 



(99) 



tan/?> 6.6 x 10" 7 , 
tanyS> 6.6 x 1(T 8 , 
tan/? > 6.6 x 10" 9 , 



if m v - leV, 
if m v - \0~ 2 eV, 



if m v - 10" 4 eV. 



(100) 



The lower limit on tan /? as a function of light neutrino mass scale m v is shown in Fig.3. 



tan/3 



6-10 



5-10" 



4 • 10 



3-10 



2-10 



1-10" 




m v (eV) 



Figure 3: The lower limit on tan/? as a function of light neutrino mass scale m v provided by 
the muonium-antimuonium oscillation experiment. The area above the curve is allowed by the 
experiment results. 



21 



Notice that the ratio of the two Higgs VEVs tan/3 are related to the light neutrino masses in the 
above inequality, although the graph (g) and (h) don't involve any neutrinos in the intermediate 
states. This results since we are using a specific model where the neutrino masses are generated 
by see-saw mechanism m v ~ 0(m 2 D IMn). The sneutrino mixing matrix is approximated in term 
of nio/MR and the heavy sneutrino masses are also of order Mr. If we take mj) to be of order 
Mw, the heavy sneutrino masses Mr in the contribution of graph (g) and (h) can be expressed in 
term of the light neutrino mass scale m v . This explains the appearance of the parameter m v in the 
inequality Eq.(99). 

However, for non-infinitesimal values of tan /3, this contribution is very small compared with 
the maximum of the second term in Eq.(83). For instance, taking the neutrino mass m v to be leV 
and assuming tan/3 > 10 -4 , the contribution of graph (g) and (h) is 

3g 4 m 2 (l + tan 2 /?) 2 

2 — 4- < 6.5 x \0~ ll GeV- 2 . (101) 

2048007r 2 M 4 v tan 4 /3 

Thus, except for the case of very small tan/3, the second term in Eq.(83) is the dominant contri- 
bution for a wide range of the parameters and its maximum is roughly two orders of magnitude 
below the sensitivity of the current experiments. 

6 The constraints from the muon and electron anomalous magnetic 
moment experiments 

One has to be careful about other constraints on the model parameters. Examples of such potential 
constraints come from the measurements of the muon and electron anomalous magnetic moments. 
The correction to the muon anomalous magnetic moment in the model under consideration is 
found by calculating the one-loop graphs shown in Fig.4. 
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(a) (b) (c) 

y j y 

Ml / \ Ml Ml / \ Ml Mr / \ Mr 

i x / s I x 

\ \ \ 

/ / / 

\ \ \ 

/ / / 

/ h \ / » \ / » \ 

^ H H MM M 

(d) (e) (/) 



Figure 4: The Feynman graphs contributing to the muon anomalous magnetic moment beyond the 
Standard Model. 

The muon anomalous magnetic moment contributed from the above graphs is 

~Mf Z (W*'c*j 

1W1 a=4,5,6 

4 

a=3 



( B5M _ 



,2 W 2 f 2 



2 2 



W~ Va=l 



+ 



(mf'mfhm 2 ( 2 4 

Z Z^> 2 ^> + Z^ c fe.> 

i=e,)i,T ' h~ \a=l a=3 

? ? ,2 2 



327r 2 M^ W 32;r 2 M 2 



=3 

2 2 2 2 2 2 

32n 2 M 2 ~ iW) 8n 2 M 2 lW ' 



B 



2 2 

mi mi 

= jjjjS % = and 

Vfi B 
F W {XyJ= f d 

Jo 



where sn - t#-, 

W° B 

C l -4x 2 (l + x) - 2x M ■ x 2 (x - 1) - 2x Va (2x - 3x 2 + x 3 ) 
Jo x M ■ x 2 + (1 - x^)x + x Va (\ - X) 



(102) 



(103) 
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F C (yv, a ) = Z-rp. ~ 4 , (104) 

6(i - yv m r 



si - 6s? + 3sr, + 6sn Insn +2 

F N (sn ) = £ . (105) 

* 6(1 -s^f 

With assumption that Mr is the largest mass scale, the dominant contribution of the graphs in 
Fig.4 to af M is 

BSM^ Z\ m l< , £ m l ( F C (V s + F C, A 

2 2 2 2 2 2 

The second and the last three terms are all decreasing functions of slepton ,chargino and neutralino 
masses. We can use the experimental bounds in Tab.3 to calculate the maximum values of these 
terms yielding 

2 2 



•(f^^ + F ^)) < 2.5x10" 



32n 2 Ml 

w- 

2 2 

mi 



g 2"> pN( \ < ,n x]0 -10 

32^/ ^ * L9Xl ° ' 

2 2 

J^—F N (t~ ) < 1 9x 10" 10 

2 2 

^Hr^(fe) * 2.1 xHT 10 . (107) 
8/r 2 M 2 ^ 



The maxima of these terms are all about one order of magnitude smaller than the present experi- 
mental bound on the contribution to a u - ^(g - 2) beyond the standard model ifToll : 

6a M = a* p -a^ = 2xl0~ 9 . (108) 

The first and third term both depend on the light neutrino mass scale m v and get suppressed. For 
instance, using the assumptions Eq.(85) and Eq.(86), the first and third terms are 

g\mlm 2 D _g\mlml ^ 26xl0 -30 



4n 2 M 2 w M 2 R \n 2 M% 
g\mlm 2 D {\ + tan 2 /?) _ gfofofll + tan 2 0) 3 3 x 1Q -3i 1+tan 2 /? 
327r 2 M 2 ,M 2 tan 2 /3 ~ 32n 2 M* v tan 2 p ~ ' X tan 2 /? 
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The first term is negligible compared with the terms in Eq.(107). However, the third term can be 
large if tan /? is very small. Therefore, the experimental bound on the muon magnetic moment will 
provide an inequality on tan f3 and m v , wihch is given by 

g 2 m 2 m 2 v (\ + tan 2 j3) 

2 M — A < 2 X 1(T 9 . (1 10) 

327r 2 M^tan 2 /3 

This inequality translates into a lower bound on tan fi as a function of the light neutrino mass scale 
m v as shown in Fig.5. 

tan/3 




m v (eV) 



Figure 5: The lower limit on tan/3 as a function of light neutrino mass scale m v provided by 
the muon anomalous magnetic moment experiment. The area above the curve is allowed by the 
experiment results. 



The electron anomalous magnetic moment beyond the Standard Model is contributed by the 
six Feynman graphs displayed in Fig. 6. 
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Figure 6: The Feynman graphs contributing to the electron anomalous magnetic moment beyond 
the Standard Model. 



The experimental bound on the contribution to the electron anomalous magnetic moment be- 
yond Standard Model is El 

6a e = a/ p -a s e M = 1.4 X 10" 11 . (Ill) 

In analogy to the muon case, this experimental limit will also generate an inequality relation of m v 
and tan yS given by 

glm 2 e m 2 v (\ +tan 2 /?) n 



32/r 2 M^tan 2 /? 
This inequality is illustrated in Fig.7. 



< 1.4x10"". (112) 
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Figure 7: The lower limit on tan/3 as a function of light neutrino mass scale m v provided by the 
electron anomalous magnetic moment experiment. The area above the curve is allowed by the 
experiment results. 

From Fig.2, Fig.5 and Fig.7 we see that for the model we are considering, the muonium- 
antimuonium oscillation experiment gives a more stringent constraint on tan j3 than the muon and 
electron anomalous magnetic moment experiments. 

7 Conclusions 

We have calculated the effective coupling constant of the muonium-antimuonium oscillation pro- 
cess in the Minimal Supersymmetric Standard Model extended by inclusion of three right-handed 
neutrino superfields where the required lepton flavor violation has its origin in the Majorana prop- 
erty of the neutrino and sneutrino mass eigenstates. For a wide range of the parameters, the 
contribution of the graphs mediated by the sneutrino and winos, W~, is dominant. The maximum 
of this contribution to the effective coupling constant is roughly two orders of magnitude below the 
sensitivity of current muonium-antimuonium oscillation experiments. Moreover, these parameters 
are not constrained by the [i — > ey decay bounds, because the lowest order contribution to the 
[i — > ey process that involves only sneutrino and winos in the intermediate states is forbidden in 
the model. Finally, there is very limited possibility that the contribution of the graphs mediated by 
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sneutrinos and Higgsino h~ is dominant if tan/3 is very small. In this case, the contributions can 
even be large enough to reach the present experimental bound. Therefore, the experimental bound 
can provide an inequality on the model parameters, which can be translated into a lower bound on 
tan /3 as a function of light neutrino mass scale m v . The constraints from the muon and electron 
anomalous magnetic moments were also investigated. For this model, the muonium-antimuonium 
oscillation experiments give the most stringent constraints on the parameters. 
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